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Summary. Many real discrete optimization problems (DOPs) are A'^P-hard and 
contain a huge number of variables and/or constraints that make the models in- 
tractable for currently available solvers. Large DOPs can be solved due to their 
special structure using decomposition approaches. An important example of decom- 
position approaches is tree decomposition with local decomposition algorithms using 
the special block matrix structure of constraints which can exploit sparsity in the 
interaction graph of a discrete optimization problem. In this paper, discrete opti- 
mization problems with a tree structural graph are solved by local decomposition 
algorithms. Local decomposition algorithms generate a family of related DO prob- 
lems which have the same structure but differ in the right-hand sides. Due to this 
fact, postoptimality techniques in DO are applied. 



1 Introduction 

Discrete optimization (DO) problems arise in various application areas such as 
planning, economical allocation, logistics, scheduling, computer aided design, 
and robotics. Application areas od discrete optimization models of OR also 
include supply chain design and management, network optimization, telecom- 
munications, VLSI routing, manufacturing, transportation, scheduling, and 
finance. The tremendous attention that DO particularly has received in the 
literature gives some indication of its importance in many research areas. 
Unfortunately, most of the interesting problems are in the complexity class 
A'^P-hard and may require searching a tree of exponential size (if P ^ NP) in 
the worst case. Many real DO problems contain a huge number of variables 
and/or constraints that make the models intractable for currently available 
DO solvers. 
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One of the promising approaches to cope with NP-hardness in solving DO 
problems is the construction of decomposition methods [31] . Decomposition 
techniques usually determine subproblems the solutions of which can be com- 
bined to create a solution of the initial DO problem. Usually, DO problems 
from applications have a special structure, and the matrices of constraints 
for large-scale problems have a lot of zero elements (sparse matrices). The 
nonzero elements of the matrix often fall into a limited number of blocks. The 
block form of many DO problems is usually caused by the weak connectedness 
of subsystems of real- world systems. 

The search for graph structures appropriate for the application of dynamic 
programming caused a series of papers dedicated to tree decomposition re- 
search (|35], [3; [28], [H], [2, [2]> [2Z])- Tree decomposition methods aim 
to merge variables such that the meta-graph is a tree of meta- nodes. Tree 
decomposition and the related notion of a treewidth (Robertson, Sey- 
mour }35| ) play a very important role in algorithms, for many A'^P-complete 
problems on graphs that are otherwise intractable become polynomial time 
solvable when these graphs have a tree decomposition with restricted maximal 
size of cliques (or have a bounded treewidth [7], [9]). 

Most of the works based on tree decomposition approach only present 
theoretical results [17], see the recent survey by Hicks et AL. But only 
few papers on applications of this powerful tool in the area of DO exist [29] . 
[24]. 

The algorithmic importance of the tree decomposition was caused by re- 
sults of Courcelle [T^ and Arnborg et al. [2] which showed that several 
A'^P-hard problems posed in monadic second-order logic can be solved in poly- 
nomial time using dynamic programming techniques on input graphs with 
bounded treewidth 

Tree decomposition based algorithms demonstrated their efficiency on solv- 
ing frequency assignment problem Koster et al. [29], ring- routing problems 
[TU] , and traveling salesman problem [TT] . 

Efficiency of tree decomposition based algorithms crucially depends on in- 
teraction graph structure of the DO problem, so that it has a time complexity 
0{n ■ 2*™+-'^), where tw is the treewidth of the graph. If the interaction graph 
is rather sparse or, in other words, it has a relatively small treewidth, then 
complexity of the tree decomposition algorithm is reasonable. 

Necessity of reduction of enumeration while solving problems correspond- 
ing to meta-nodes of the tree decomposition causes expediency and an urgency 
of development of tools that could help to cope with this difhculty 

In this paper, discrete optimization problems with a tree structural graph 
are solved by local decomposition algorithms that belong to dynamic program- 
ming paradigm. Local decomposition algorithm generates a family of related 
DO problems which have the same structure but differ in the right-hand sides. 
Due to this fact, postoptimality techniques in DO are applied. 
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2 Discrete optimization problems with constraints and 
their graph representations 

Consider a DOP with constraints: 

nmx/(X) = max^/,(X'^), (1) 

subject to 

A^s.Xs, < h, ieM = {1,2, . . . ,m}, (2) 
=0,1, je7V = {l,...,n}, (3) 

where 

X = {xi, . . . ,Xn} is a set of discrete variables, functions fiiX"^) are called 
components of the objective function and can be defined in tabular form, 
X'' C X, k £ K ~ {1,2, . . . ,t} , t is the number of components of objective 
function, is a set of indices of components; 

5,C{l,2,...,n}, zeM. (4) 

We shall consider further a linear objective function ([5]): 

n 

F{xi, . . . ,Xn) — F{X) ~ CnXn ~ "^^CjXj max (5) 

i=i 

Definition 1. f^. Variables x £ X and y £ Y interact in DOP with con- 
straints if they both appear either in the same component of objective function, 
or in the same constraint ( in other words, if variables are both either in a set 
X^ , or in a set XsJ- 

Graph representation of a DOP structure may be done with various de- 
tailization. Structural graph of a DOP defines which variables are in which 
constraints. 

An interaction graph [5] represents a structure of the DOP in a natural 
way. 

Definition 2. J^]. The interaction graph of a DOP is an undirected graph 
G — {X, E), such that 

1. Vertices X of G correspond to variables of DOP; 

2. Two vertices of G are adjacent iff corresponding variables interact. 

Further, we shall use the notion of vertices that correspond one-to-one to 
variables. 

Definition 3. The set of variables interacting with a variable x € X , is de- 
noted by Nb{x) and called neighborhood of the variable x. For corresponding 
vertices a neighborhood of a vertex v is a set of vertices of interaction graph 
that are linked by edges with v. Denote the latter neighborhood as Nbciv). 
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Let 5 be a vertex set of the graph. Introduce the following notions: 

1. Neighborhood of a set 5 C V, Nb{S) = U^gs ^b{v) - S. 

2. Closed neighborhood of a set S" C y, Nb[S] = Nb{S) U S. 

3. If S* = {ji, . . .,jg} then Xs = {xj,,. . .,XjJ. 



Example 1. 



2xi + 3x2 + Xs + 5x4 + 4a;5 + Gxq + xr 
3xi + 4X2 + X3 

2x2 + 3^3 + 3x4 

2x2 + 3x5 

2x3 + 3a;6 + 2x7 

Xj = 0, 1, j = 1, . . . , 7. 



■+ max 

< 6, (Ci) 

< 5, {C2) 

< 4, (C3) 

< 5, (C4) 




Fig. 1. Interaction graph for example. 



We need following notions. A clique is a set of vertices that induce a 
complete subgraph of G, and a maximal clique is a clique which is not a 
subset of any other clique. A clique C in graph G is maximal, if G is not a 
subset of any other clique in G. A spanning tree of a graph is a tree that 
includes every vertex in the graph. The clique graph of G, is the intersection 
graph of the family of maximal cliques of G. 
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2.1 Tree decomposition in discrete optimization. Chordal graphs. 
Triangulation 

Dynamic programming is a very powerful algorithmic framework in which 
an optimization problem is solved by identifying a set of subproblems and 
solving them one-by-one, smallest first, storing solutions and using the stored 
solutions to small problems to find solutions of larger ones (recursively defining 
value of optimal solution), until the whole set of them are solved. Dynamic 
programming computes recurrences efficiently by storing partial results in 
tables. 

A tree is a recursive data structure because each child of a node in the tree 
has a tree of its descendants. Due to this fact, many of the important algo- 
rithms to access and manipulate trees are easily expressed using recursion and 
particularly dynamic programming which computes recurrences efficiently by 
storing partial results in tables and hence can be effectively used for solving 
optimization problems on trees. Dynamic programming works best on objects 
which are linearly ordered, e.g. the left-to-right order of leaves in a tree. Dy- 
namic programming starts at the leaves of the tree and proceeds from smaller 
to larger subproblems (corresponding to subtrees) that is to say, bottom-up 
in the rooted tree. 

One reason why many optimization problems that are hard on general 
graphs are easy on trees is that trees do not contain cycles. 

Above mentioned facts and an observation that many optimization prob- 
lems which are hard to solve on general graphs are easy on trees makes the 
detection tree structures in a graph a very promising task. A powerful tool of 
algorithmic graph theory like tree decomposition f35] can help detect trees 
and obtain the treewidth, a measure of the "tree-likeness" of the graph. 

Using dynamic programming techniques on tree decompositions of bounded 
treewidth, hard optimization graph problems can often be solved in polyno- 
mial time [3 . 

The notions of treewidth and tree decomposition were introduced by 
Robertson & Seymour in their seminal paper 155] on graph minors. The 
best known complexity bounds are given by the treewidth tw (Robertson, 
Seymour [3S]) of an interaction graph associated with a DOP. This parame- 
ter is related to some topological properties of the interaction graph. It leads 
to a time complexity of 0{n ■ 2*"'+i). Tree decomposition based methods aim 
to merge variables so that the resulting meta-graph is a tree of meta-nodes. 
A more detailed introduction to tree decompositions is given in [28 and in 
surveys |S] , [23] , [JDj . Most of the works based on tree decomposition approach 
only present theoretical results [17]. 

Definition 4. Let G — {V, E) he a graph. A tree decomposition of G is a 

pair (T; X) with T = (/; F) a tree and X — {Xi | / G /} a family of subsets 
of V , one for each node of T , such that 
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(a) for every edge (m, v) there is an i ^ I with u ^ Xi, v £ X^, 
(iii^running intersection property) for alii, j, I G /, ifi < j < I, then XiHXi C 
X,. 

Remark 1. To construct a tree-decomposition we merge the vertices of G to- 
gether to form meta-nodes ("supernodes" [1] or "bags" [H]); each meta-node 
is a subset of the vertices of G and we connect these meta-nodes to form a 
rooted tree T. The meta-nodes do not have to be disjoint, i.e., there might be 
nodes of the graph contained in more than one meta-node. This grouping and 
connecting has to be done in such a way that - for all edges e of the graph: 
there exists a meta-node containing both cndpoints of edge e, and - for all 
vertices v of the graph: all meta-nodes containing v together with the edges 
between those meta-nodes in T form a connected subtree of T. 

The width of a tree-decomposition is the number of graph- vertices of the 
largest meta-node minus 1. 

As finding an optimal tree-decomposition is A'^P-hard, approximate optimal 
tree decompositions using triangulation of a given graph are often exploited. 
Given a triangulated (or chordal) graph, the set of its maximal cliques cor- 
responds to the family of subsets associated with a tree-decomposition (so 
called clique tree [6]). When a tree-decomposition is exploited, usually one 
considers approximations of optimal triangulations by clique trees [26]. Hence, 
the time complexity is then 0(n -T" with w+ + 1 the size of the largest 
cluster (ui -|- 1 < -I- 1 < n). The space complexity is 0(n • s • 2**) with s the 
size of the largest minimal separator |27j . 

Definition 5. A clique tree of G is a tree T = {K, E) whose vertex set is 
the set of maximal cliques of G such that each of the induced subgraphs T[Kv\ 
is connected. 

Definition 6. A graph is chordal (triangulated, perfect elimination graph, 
rigid circuit 1151, monotone transitive 1361) if every cycle of length > 3 has a 
chord (i.e., an edge joining two nonconsecutive vertices of a cycle). 

All induced subgraphs of a chordal graph are also chordal. If G is a chordal 
graph, then any clique tree of G is also a tree decomposition of G. However, 
the converse is not necessarily true. 

Theorem 1. [1 7/ Let G = (V, E) he an undirected graph, and let K he the set 
of maximal cliques of G, with Ky the set of all maximal cliques that contain 
vertex V ofG. The following statements are equivalent: 

(i) G is chordal. 

(ii) G is the intersection graph of a family of subtrees of a tree. 

(iiijThere exists a tree T — (K,E) whose vertex set is the set of maximal 
cliques of G such that each of the induced subgraphs T[Ky] is connected. 
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Clique graphs are not very useful on general graphs, since these can contain 
n\/{kl{n — k)l) different cliques of size k. It follows from this that clique graphs 
of general graphs can be exponentially large, and it is no surprise that finding 
the maximal clique for general graphs is hard. Chordal graphs on the other 
hand have limitations that make clique graphs useful. 

Lemma 1. (T)IRAC JlSf ) A chordal graph G contains at most n maximal 
cliques. 

Gavril TTj proved that every chordal graph can be represented by a clique 
tree limiting the number of edges to 0{n). Indeed, from theorem [T] it follows 
that for the chordal graph G there exists a tree T = {K, E) whose vertex set 
is the set of maximal cliques of G such that each of the induced subgraphs 
T\Ky\ is connected. From Lemma [T] and the fact that T is a tree, follows that 
a clique tree has at most n nodes and n — \ edges. 

Theorem 2. (13ernstein and Goodman maximum weight span- 

ning tree of the clique graph of a chordal graph G is a clique tree of G. 

Theorem 3. (Ho AND Lee \25^ ) Given a chordal graph G and a clique tree 
T of G, a set of vertices S is a minimal separator of G iff S ~ Gi H Cj for an 
edge (Gi, Gj) in T . 

Corollary 1. fHo and Lee J2^ ) A chordal graph G has at most n — \ min- 
imal separators. 

These results show how to build tree decompositions using elimination game 
algorithm which triangulates an initial interaction graph. For triangulated 
graphs it is rather simply to find maximal cliques and to build the clique tree. 

2.2 Elimination game and tree decomposition 

The process of interaction graph transformation known as Elimination Game 
was first introduced by Farter [33] as a graph analogy of Gaussian elimina- 
tion. The input of the elimination game is a graph G and an ordering a of G 
(i.e. a(v) = i li V \s z-th vertex in ordering a). Elimination Game according 
to [23] consists in the following. At each step i, the neighborhood of vertex 
Xi is turned into a clique, and Xi is deleted from the graph. This is referred 
to as eliminating vertex Xi. The filled graph G+ = {V,E^) is obtained by 
adding to G all the edges added by the algorithm. This resulting graph G+ is 
a triangulation of G (FULKERSON & GROSS [16]), i.e., a chordal graph. 

Different filled graphs result from processing the vertices of G in different 
orders. Thus in order to find a low fill-in, it is necessary to find a good order 
on the vertices of the given graph before running elimination game. Finding 
an ordering that results in the minimum fill-in is an A^P-complete problem 

Ha- 

Elimination Game can also be implemented so that a is generated during the 
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course of the algorithm. In this case, we can at each step i choose a vertex 
V of the ehmination graph G*~^ according to any desired criteria, and set 
a{v) = i, to define an ehmination ordering a. One well known heuristic called 
Minimum Degree chooses a vertex v of minimum degree in G"~^ at each step 
i. 




a) Initial interaction graph; b) after elimination . 




c) after elimination . d) after elimination x^ and x^ . 





e) after elimination x^ and x, . 

Fig. 2. Elimination Game. 

The procedure to solve an optimization problem with bounded treewidth 
involves two steps: 

(i) computation of a (good) tree decomposition, and 

(ii) application of an (dynamic programming) algorithm that solves instances 
of bounded treewidth in polynomial time. 

To describe how tree-decompositions are used to solve problems with dynamic 
programming, let us assume we find a tree-decomposition of a graph G |39j . 
Since this tree-decomposition is represented as a rooted tree T, the ancestor/ 
descendant relation is well-defined. We can associate to each meta-node X the 
subgraph of G made up by the vertices in X and all its descendant meta-nodes, 
and all the edges between those vertices. Starting at the leaves of the tree T, we 
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can compute information typically stored in a table, in a bottom-up manner 
for each bag until we reach the root. This information is sufficient to solve 
the subproblem for the corresponding subgraph. To compute the table for a 
meta-node of the tree-decomposition, we only need the information stored in 
the tables of the children (i.e. direct descendants) of this node. The problem 
for the entire graph can then be solved with the information stored in the 
table of the root of T. 



3 Local decomposition algorithm for discrete 
optimization 

During the study of complex objects it is not always possible (and expedient) 
to obtain (or to calculate) complete information about the object as a whole; 
therefore is of interest to obtain information about the object, examining it 
in parts, i.e., locally. 

Yu.I. Zhuravlev [44 introduced and investigated local algorithms for 
calculating the information about properties of objects. The local algorithm 
can be described as follows. For a given set of sets {m} , m = {Ui} , i = 
1, 2, . . . , |m|, and for each element U d m let us determine a neighborhood 
S{U, m) in to; these neighborhoods should satisfy the following conditions: 

• U<ES{U,m); 

• S(U, to) C to; 

• If [/ G TOi, U e TO2, S{U,m) □ TO2 C TOi then S{U,mi) — 5(L/, TO2). 

Algorithm A is completely determined by the set of predicates Pi , . . . , P; , 
by the partition of this set into a subset of basic predicates Pi , . . . , Pr 
and auxiliary predicates Pr+i,...,Pi , by the set of monotonic functions 
(fii, . . . ,ipi, ipi ~ (pi{U, ai, . . . ,ai, S, TO*) and by the ordering algorithm A-^. 

LA uses a dynamic programming paradigm and computes optimal partial 
solutions of subproblems that correspond to the blocks of the DOP. 

Each step of the LA 21 38j consists of changing neighborhoods and replac- 
ing index p with p+1 (although it is possible to pass, also, from Sp to Sp+p); 
for each fixed assignment of the variables of boundary ring the values of the 
variables of the corresponding neighborhood are stored, in this consists one of 
the important differences between LA 21 and LA A: an information not only 
about the predicates, but also about the values of variables is memorized; this 
information is called as an indicator information. 

Let us consider the LA 21st for the solution BT ILP problems (O, 
([3]), where the matrix A has BT structure with the tree D which contains of 
k blocks. 

Consider a vertex r of the tree D and introduce a tree Dr which consists 
of the vertex r and all its descendants. 
Introduce the necessary notations: 
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• S'r is a set of indices of variables which belong to block B^; 

• S^^' is a set of indices of variables which belong simultaneously to blocks 
Br and B^' ; 

• if S = {ji, . . .,jq}, then Xs = {xj,,. . . 

• Pr is a vertex-ancestor for the node r; 

• Jr is a set of descendants of the node r. 

Consider a LA [38^ for solving DO problems with a tree structural graph, i.e., 
problems in which it is possible to find the set of the neighborhoods of different 
variables so that one variable can belong to two neighborhoods only and the 
graph of intersections of these neighborhoods is a tree. The LA solves this 
DO problem, moving bottom-up, i.e., from the neighborhoods corresponding 
to leaves of the tree, to the neighborhood corresponding to the root of the 
tree T. Let Bi = Ui) , B2 = {S2, C/2) , ■ ■ • , Bk = {Sk, Uk) be a set of 
the blocks (neighborhoods) of some indices ji, ■ ■ ■ ,jk of some variables, where 
Sr , Ur are, respectively, the sets of the indices of variables and constraints for 
the rth neighborhood, r = 1, . . . , fc and 

k 

[jUr^M^{l,...,m}, (6) 

fc 

[jSr^N^{l,...,n}, (7) 

Ur.nUr, =0, ri ^ra, (8) 

Sri ^ Sr2 n S'r., = for any triple of different indices ri, r2, ra. (9) 

Consider a node r of the tree T and define a tree containing the node r 
and its children. 

Introduce notions: 

• Sr is a. set of indices of variables that are in the block Br\ 

• Srr' is a set of indices of variables that are in the blocks Br and B^i , i.e., 

Srr' — Sr P| Sr' ', 

• Pr is a, node-parent of the node r; 

• is a set of children of the node r. 

Then Xs is a meta-variable consisting of variables common for blocks Bp^^ 
and Br (here Sp^r = Sp^ P| Sr)- 

Denote as Zt^ the following problem: for each assignment Xs^^^ to find 
^Sr- > , such that 

hsAXs^J - fo^ (^s„.0 = ^ max J Cs^Xs^ + V IfoA^s ,) + Cg , Xs 
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y r' 

subject to 

As^Xs^<K~ J2 M^r'^s^.' -As.^rXs,^.- 

r' GJr 

Here fT^{Xs^^^) is an objective value of subproblem corresponding to the 
tree T^, fr , [Xs , ^ — Ct , Xt , . It is possible to assign this value to the 

root B^i of the tree Tr and write: hs , {X^^ fl s ' ) ^ , (^Xs , ^ . 

It is easy to see that if we fix a vector Xs ^ , then the problem is decom- 
posed into two problems: the first one corresponds to the tree T^; and the 
second one to T — Tr. An application of LA 21bt for solving DP problems 
with a tree structural graph is based on this property. 



4 Solving of a concrete DOP with the finding of tree hke 
structure and applying of local decomposition algorithm 

Consider the DOP of example 1. 

Finding of the tree structural graph (tree decomposition) 

In Fig. [2] results of elimination game algorithm are shown. Since during 
elimination process new fill-in edges are not added, the elimination game pro- 
cess is equivalent to searching simplicial vertices and corresponding maximal 
cliques. In Fig. [3] these maximal cliques and links between them are shown. 
Local decomposition algorithm can be applied to this clique tree. Other pos- 
sible way of finding of the clique tree is using of maximal spanning tree in the 
dual graph. 

Applying the local decomposition algorithm to DO problem 

Let us solve the subproblem corresponding to the block Bi . Since this block 
is adjacent to the block B4, we have to solve DOP with variables Xbi-b^ for 
all possible assignments f| 54. Thus, since Xbj_s4 = {2:5} and p| = 
{0:2}, then induced subproblem has a form: 

hBi{x2) = max {4x5} 

subject to 

2x2 + 3a;5 < 4, Xj = 0, 1, j £ {2, 5} 



Solution of the problem can be written in tabular form: 
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Fig. 3. Tree decomposition for tlie example 1. 

Table 1. Calculation of hBi{x2) 



X2 




XI{X2) 





4 


1 


1 









Next subproblem corresponding to a leaf (or meta-node) B2 of the clique 

tree is 
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Hb^ {x^) = max {Gxg + 2:7} 

S.t. 

2xs + 3a;6 + 2x7 < 5, Xj = 0, 1, j e {3, 6, 7} 
Solution of this subproblem: 

Table 2. Calculation of /i 32(3^3) 



X3 


hB2 


XI{X3) X^ix^) 


1 


1 


7 1 


1 





6 1 



Subproblem corresponding to the block S3 has the form: 
hB3{x2,X3) = max {/iB^ (2:3) + ^xa} 

X4 

S.t. 

2x2 + 3x3 + Sxi < 5, Xj = 0, 1, j e {2, 3, 4} 
Table 3. Calculation of hB3{x2,X3) 



X2 X3 




XI{X2,X3) 





12 


1 


1 


6 





1 


12 


1 


1 1 


6 






The last problem left to be solved is: 

hBi= max {hBi{x2) + hB3{x2,X3) +2xi} 

Xl,X2,X3 

S.t. 

3x1 + 4a;2 + ^3 < 6, Xj = 0, 1, j G {1, 2, 3} 
Table 4. Calculation of /1B4 





Xl X2 X3 


18 


1 



The maximal objective value is 18. To find the optimal values of the vari- 
ables, it is necessary to do backward step of the dynamic programming pro- 
cedure: from table 4 we have x^ = 1, X2 = 0, = 0. From table 3 using the 
information X2 = 0, x| = we find x| = 1. Considering table 2 we have for 
xl—0: xl = l,Xj = 1. From table 1 we find for x^ = 0: X5 = 1. The solution 
is (1, 0, 0, 1, 1, 1, 1); maximal objective value is 18. 
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5 Postoptimal analysis and local algorithms 

5.1 Postoptimal analysis in DO 

Decomposition and sensitivity analysis in DO are closely related. Sensitivity 
analysis follows naturally from the duality theory. Decomposition methods 
consist of generating and solving families of related DO problems that have 
the same structure but differ as the values of coefficients. Sensitivity analysis 
allows using information obtained during solving one DO problem of the fam- 
ily of related DO problems in solving other problems of this family. Due to the 
lack of full-fledged duality theory in DO, sensitivity analysis for DO problems 
is not sufficiently developed [19j, [31]. A number of useful tools of sensitivity 
analysis in DO are derived for integer programming in [19j . A technique of 
sensitivity analysis proposed in '37] computes a piecewise linear value function 
that provides a lower bound on the optimal value that results from changing 
the right-hand sides of constraints. 

Recently, an interesting application of binary decision diagrams (BDD) 
(introduced earlier in computer science community) was proposed by Hadzic 

6 Hooker [21] for the purposes of postoptimal analysis in DO. 

5.2 Postoptimal analysis in LA 

LA systematically proceeds with so called parametric DO problems [5]: 
Definition 7. A parametric DO problem is 

mm{f{X)\X e Ci, i = 1,. . . ,m; xj & {0, 1}, j = 1, . . . ,n}, 
X — p 

where Ci is a set of feasible solutions of the constraint i, z = 1,...,to. 

Thus, an optimization problem is in parametric form when the objective func- 
tion is optimized not over the entire set X, but only over a subset X — P, for 
all possible assignments of the variables of P. 

Below we show that these parametric DOPs generated by LA lead to 
the possibility of exploiting postoptimality and sensitivity tools in the LA 
procedure. Consider DOP ([S]), ([2]), ([3]) with tree structural graph and use LA 
for its solving. Then for the block Br we have to solve a family of DOPs: 

Cs^Xs^ + [Cs,_i,,Xs_, „ + C^^_,X„^_,] max 

s.t. 

These DOPs should be solved for all binary assignments Xs^_i ^ and Xs^ ^^i- 
Denote 
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It is clear, that it is better for each Xs^_^ ^ and Xs^^^j^^ to solve a problem: 

Cs^Xs, -> max (10) 

s.t. 

^ aijXj < hi (Xs,_i,jXs^,,+i) , i&Ur, (11) 

Xj = 0,1, j G Sr. (12) 

It is possible to use information obtained during the solution of some DO 
problems of the family (fl0l[T2|) for solving other problems of this family using 
postoptimality analysis (PA) [19], [20], [21], [37]. The more efficient the PA 
procedure is, the better LA will work. 

5.3 Postoptimal analysis for an implicit enumeration algorithm 

We show how to use PA in LA using as an illustrative example, the simplest 
case of an implicit enumeration algorithm (18j that generates partial solutions 
and tries to fathom them using 3 tests. 

The DO problem (Il0l)-(|l2]) with vector right hand side h {Xs^_^^^ 
will be denoted as Zr (Xs^ ^ ^\Xs^^^^\ Introduce a partial order over a set 

of DO problems {Zr (Xs,._, jAs,.^^ J }: DOP < = Zr {x's^_^ ^\X'g^ pre- 
cedes DOP Z'; = Zr (^S._,J^S.,.+ ,) if (^5._,,J^5.,.+,) < h (^l-,,J^l,.+ > * e 

Ur or graphically 





X's,. ) 


1 







Fig. 4. Partial order of DOPs. 

Let US solve the problem Z^ by implicit enumeration. Let P be a partial 
solution (PS) fathomed by the tests of the implicit enumeration. There are 3 
possible cases: 
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a) PS P is fathomed by test 2, i.e., the best completion of P is feasible. Then 
this completion is feasible for the problem too, so as 




Thus, each incumbent of Z^, is a feasible solution of Z^ ,i.e., the objective 
function value of Z^ is higher (better) than the corresponding objective 
function value of Z^. 

b) PS P is fathomed by test 1 in the problem Z^, i.e., z < z*. Then P is 
fathomed in Z^, , too, as^ < z * < z *. 

c) PS P is fathomed by test 3 in Z^. 

It is clear, that in cases a) and b) it is senseless to fathom PS P in the 
problem Z^ , because P is automatically fathomed. Thus, while fathoming PS 
in Zj. , it is interesting to study only those PS P fathomed by test 3 in Z^. 

Let us consider a family of 2\s^-i,r\+\s^,^+i\ dqPs (UHll-IISl) with nonneg- 
ative coefficients (multidimensional knapsack problems). Using partial order 
of DOPs described above it is possible to order members of this family as it 
shown in Fig. [5l Consider a process of solving the family of DOPs {Zr}. Let 
P be a PS fathomed in {Z^(ll|l)} (level in Fig. E]). If PS P is fathomed 
by tests 1, 2, we can exclude this PS P from consideration (and fathoming) 
in other DOPs. If PS P was fathomed by the test 3 in {Zj.(ll|l)}, then we 




Fig. 5. Partial order of DOPs (multidimensional knapsack). 
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pass to a DOP from the level 1, say, {Zr(01|l)} and consider the PS P in this 
problem. 

If the PS P is not fathomed by the test 3 in {Zr(01|l)}, then there exists 
one of the following cases: 

• test 1 is true, then backtrack to {Zr(ll|l)} and pass to the next problem 
of level 1; 

• test 2 is true, then also backtrack to and pass to the next 
problem of level 1; 

• both tests 1 and 2 are not true, then extend PS P : P' = {P,ji) and try 

to fathom P' . 

If any of the tests 1, 2 are true, then do usual backtracking, i.e., PS P" = 
(P, — ji) is considered. If some extension of P : P'" = {P,ji,—j2,...,jf) 

is fathomed by the test 3, then wc go to one of problem's ancestors (say, 
{.Zr(00|l)}) and try to fathom by the tests of implicit enumeration. 



Conclusion 

Local decomposition algorithms combined with tree decomposition methods 
are a promising approach that enables solving sparse discrete optimization 
problems from applications. The performance of these algorithms can be im- 
proved with the aid of postoptimality analysis. 

A promising direction of future research is the development for ef- 
ficient schemes of postoptimality analysis embedded in local decomposition 
algorithms combined with tree decomposition methods. 
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